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for the degree 
of 

MASTER OF TECHNOLOGY 
by 

B. VISHWANATH BALIGA 
Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 
AUGUST, 1986 


The present work is an attempt to analyse fracture in 
a welded asperity junction. The welded asperity junction is 
formed during rubbing motion between two surfaces, wherein 
the hi^ points of the surfaces come in contact with each 
other and due to the intense local stresses a strong welded 
junction of these asperities (high points) is formed. Adhesive 
wear, has been found to occur when these junction formed 
fracture. Wnile reviewing the available literature it was 
found that no exact analytical work on the formation of a 
wear particle as a result of fracture of asperity junction 
has been done. A few experimental results are available wherein 
the pa-'oh of frqcture of these junctions have been analysed; 

The finite element method was applied to this present 
analysis becayse it is a versatile method and it can easily 
accommodate material non-linearity, geometric non-linearity 
and complex boundary ccanditions when compared v'ith numerical 



methods^ Based on the displacement formulation, an elastic 
analysis computer code using FEM is developed to predict 
the path of fracture in a brittle junction. Also, a ductile 
junction was analysed for the fracture path using a elasto— 
plastic analysis computer code based on the initial stress 
method; Material was considered to obey Von— Mises yield 
criterion. Six node isoparametric elements are used in 
the discretization of the continuum; 

The junction is analysed at a stage in its formation, 
where it is perfectly symmetric and the fracture is just about 
to sikoixt. A fracture path is traced with the fracture occurring 
in only one of the asperities as happens in a bimetallic 
junction .fracture path is also traced with the fracture occurring 
in both the asperities simultaneously as is the likely case 
in junctions of similar metals. Three different finite element 
meshes are used to reach the path of fracture; 



CHAPTER 1 


INTRODUCTION TO THE ADHESION WEAR 
1 ,1 Introduction : 

The study of wear due to the fracture of minute asperiti 
resulting from large plastic deformations of welded asperity 
junctions formed on the sliding surfaces has been investigated 
by several research workers in this area. These welded 
asperity junctions are formed between the two metals wherever 
asperities engage and these junctions adhere strongly under 
intense local stresses. Each one of these welded asperity 
junctions goes through a life cycle; formation, pj.astic defor- 
mation and eventually fracture; new junctions being formed as 
fast as the old ones are fractured; 

The problem of fracture is the central problem in 
the Science of resistance of materials. At present when 
speaking of fracture one usually implies a study of the 
conditions under which there is a crack or a system of cracks 
propagating in a body. But cracks may be very different in 
nature and they are considered on different scales; On 
one hand, the fracture of a crystal grain begins with the 
formation of a submicroscopic crack with the separation of 
two atomic layers by such a distance that the force of inter- 
action between the atoms are zero. On the other hand we have 
the macroscopic crack (as in the case of welded turbine 
rotor or in a nuclear reactor.) where the length and width of 
the crack may be estimated by the milimetre; 
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The fundamental problem on which most investigators 
have concentrated in recent years has been the problem of 
finding the conditions for the propagation of a. crack. The 
propagation of the crack begins, when the plastic strain 
near its tip becomes large. The criterion for the initiation 
of crack propagation, which forms the basis of fracture 
mechanics does not follow from the equilibrium equations, 
but is an additional boundary condition in the solution of 
the problem,- Fracture condition is said to be reached if a 
crack like cut can propagate [1]. 

The study of crack propagation has been done by 
several research workers using numerical techniques. At 
present one of the most generally used method, having the 
least limitations is the finite element method. In the last 
few years the finite element method has found extensive 
application in the solution of various complex problems in 
structural mechanics, elasticity and other fields. In many 
situations an adequate engineering model of the system is 
obtained by replacing the system with a finite number of 
well defined components called elements. Even if the number 
of elements is very large the problem can now be solved 
readily with the advent of computers. The various discre- 
tization methods suggested from time to time to solve the 
continuum problem involve approximations which approach the 
true continuum solution as the number of discrete variables 
increase ; 
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Finite element method can be interpreted as a piece- 
wise application of the variational methods, in which the 
approximation functions are algebraic polynomials and the 
undeterfflined parameters represent the values of the solution 
of a finite number of preselected points called nodes; on 
the boundary and in the interior of the element: Finite 
element method can thus be systematically programmed to 
accommodate complex, non-linear stress - strain behaviour 
and complex boundary conditions which will be very difficult to 
accommodate in other numerical methods like the method of 
least squares, Ritz method etci 

1 «2 Review of Previous Work ; 

( a ) On Asperity Junction : 

It is a well established fact that when two solid 
bodies, in contact, are in the relative sliding motion, 
surface deterioration invariably occurs. The surface deterior- 
ation ultimately leads to wear of these bodies. Because of 
the waviness and roughness of surfaces, contact between two 
sliding bodies (Figi 1,1) occurs only at certain high points 
known as asperities. These asperities are subjected to very 
high stresses even at small loads due to their small dimensions^ 
Due to the high stresses contact of the two bodies at the 
asperities generally results in the formation of welded 
asperity j'unctions, which have been the subject of investi- 
gation for many research workers. Various experimental 
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(a) 

PARALLEL TO THE SURFACES 


FIGUDIAGRAM OFA JUNCTION® DURING THE INITIATION OF SLIDI 
(.b) DURING STEADY SLIDING 



5 


techniques like electrical resistance measurements [8,9], 
optical [10] have bean carried out, based on these it is 
generally accepted that the junctions formed by the asperity 
interactions are roughly of the order of 10 to 20M- ; 

Rabnowicz [11] used an auto^correlation analysis to 
determine the average junction size for sliding metals. 

The size distribution of the junctions is given in Fig,1.2, 

Green [12] has modelled as asperity junction (Fig,1,3)« 
At any ins^tant when the junction is formed, force acting 
through the junction has a normal component (compressive) 
and a tangential force F^ (shear) as shown in Fig; 1:3. 

Brockley and Fleming [13] have conducted experimental work 
on a model junction simulating metallic wear; They have 
predicted the wear rate by dividing the volume of the detached 
particle formed to the distance traversed during a particular 
deformation; Fig; 1;4 shows how the fracture starts during 
the process of sliding and plastic deformation; 

Literature survey conducted by the author did not 
reveal any Finite Element application to wear analysis. As 
already mentioned, the formation of a wear particle during 
adhesive wear is a complex phenomena. No exact analytical 
work to estimate the amount of material that can get lost 
by a solid during the interaction with another body has been 
done. However, a large number of empirical equations have 
been given. Only Gxem [14] and Brockley [13] have discussed 
the formation of a wear particle by tracing the fracture path 
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FIG 1.2 FRACTION OF CONTACT AREA AS A 
function ofjunction width 




FIG 1.3 MODEL OF AN ASPERITY JUNCTION 
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at a interface. Green discusses it only for a perfectly 
plastic material, Brockley through experiments has shown 
how a wear particle can be formedi 

( b) On Fracture s 

The importance of fracture mechanics for the study 
of crack propagation and its effects in structures has 
increased considerably over the recent years. A factor 
known as the stress intensity factor, K, has been developed 
as the critical value for crack to. propagate. The magnitude 
of K depends op the structural geometry and the loading 
system. If the value of K exceeds a certain critical valuel, 

K ., crack instability is assumed to exist, so that K is 
now effectively a material property which represents the 
resistivity of the material to fracture [ 1 ] ; 

A number of solutions for K, have been deduced for 
various crack sizes and loading for relatively simple shaped 
structures, Paris and Sih [2^ give a comprehensive handbook 
of such results. However, for more realistic complex shapes 
encountered in practice more general techniques are required. 
Use of Finite Element Method enables such computations to be 
made, Jerranj and Hellen [3], using finite element have derived 
stresses and displacements around the crack tip and then 
using the computed values have evaluated K, The computed 
values of stresses and displacements substituted into the 
known crack tip field equations which relate stress and 
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displacement to K, yield the value for the stress intensity 
factor; 

Watwood[4] has done considerable work based on the 
energy release criteria of Griffiths [ 5] using this method 
Watwood and others have calculated the energy differences 
between successive points along a mesh representing a growing 
crack and used this to find the possible path of fracture. 

This method is sometimes referred to as virtual crack 
extension method [ 6 ]i 

A different method has been developed by Parks [7 ]; 

He uses the J-integral approach on a piecewise straight path 
around the crack tip along the sides of the element used. 

Parks [7] has also shown that using a constant stress .triangleai 
the virtual crack extension method is equivalent to the 
J-integral approach when expressed in terms of stresses and 
displacements; 

1 .3 Junction Formation between Contacting Bodies ; 

No matter how plane the surface of a solid may 
appear to be, when examined at microscopic level it is always 
wavy and rough; The lengtji of the waves varies from 
1 mm - 10 mm and their heights vary frcm 20^50P'[143 ; A 
profile of a actual surface is shown in Fig;1.5 [15]; 

Due to this distortions in shape, the actual surfaces are 
in a position to develop contact only over small discrete 
areas; 
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FIG. 1.S ROUG HNESS AND WAVINESS OF SURFACES 

(A) PRCFIlf OF A TURNED COPPER SURFWEClOltpi STORIED MV&NinC4T]ONCBELOW)E6lUAL MA&NIfICAT] 

(B.^ WAV1NE5S DIAGRAM OF A TURNED SURFACE 
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A closer examination of these discrete regions, of 
contact has revealed further that the actual contact between 
bodies occurs at certain high points located within these 
regions, these regions are called as contact points or more 
commonly asperities. Problems connected with metal transfer 
and wear during sliding process are very much dependent on 
the nature of phenomenae taking place at these contact asperi- 
ties, These asperities are subjected to extremely high stresses 
even at ordinary loads, because of their very small dimensions; 
Due to the high stresses contact of the two bodies at the asperi- 
ties generally results in formation of junctions which have been 
the subject of investigations for many research workers. Various 
experimental techniques like electrical resistance measurement 
[16, 17], optical [18] and electron microscopy [19] have been 
employed to study these asperity junctions; 

Rabinowitz proposes the size distribution of the 
junction to be given as in Fig;1 ;2, Ling [20], Greenwood and 
Wiliams on [21] and many others have also analysed surface 
profiles using statistical approaches and ha ve^ generally 
concluded that the junction size varies from 5-20 P under 
normal conditions of rubbing; 

Adhesive Wear ; 

The tendency of the contacting surfaces to adhere 
arises from the attractive forces which exist between the 
surface atoms of the two materials. If two surfaces are 
brought together and then separated, either normally or 



13 


tangentially^ the attractive forces act in such a way as to 
attempt to pull material from one surface on to the other. 
Whenever, material is removed from its original surface in 
this way an adhesive wear fragment is created. The strength 
of bonding at the points of adhesion is sometimes so great 
that while attempting to free the surfaces separation takes 
place not along the interface but in one of the bodies 
itself resulting in metal transfer and subsequent metal 
removal, this will occur if the force required to break 
through the interface of the materials is larger than the 
force required to break through some continuous surface 
inside one of the materials. It has been experimentally 
observed that because of the intense local stresses at the 
junction, strong adhesion generally occurs at the interface 
and hence fracture takes place mostly along some continuous 
path inside one of the materials. Thus under some conditions, 
loss of material may occur from one or from both of the 
rubbing pair. This type of wear is called as adhesive wear. 
This is the most conmon form of metal wear and exists predo- 
minantly at moderate cutting conditions. The removal of 
material takes place in the form of small particles which 
are generally transferred to the other surface but may also 
come of in loose form; 

The way in which wear particles form is graphically 
illustrated by the experiments of Green [12] who used two- 
dimpnsional models of various metals and plasticene to denote 
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asperit:.es, and then sheared the asperities; 

The stages through which junctions have been observed 
to be passing can be understood with the aid of Fig. 1 .6; 

The figure illustrates the deformation of asperities I and 
II of similar hardness. After the initial adhesion mgst of 
the centre of the junction is sheared uniformly? Fig.1;6(c), 
Due to this shearing the interface gets stretched. This 
would help to break any surface films present and thereby 
strengthen adhesion. As the tendency for the asperities to 
pass each other continues an intermediate stage shown in 
Fig;1 ,6( d) is reached when the junction has acquired a 
symmetric shape. Since the adhesion is assumed to be 
fairly strong the junction deforms as a single body and thus 
necking and fracture occurs in the most critically stressed 
region of the junction. In this manner a small particle is 
transferred from one asperity to the other one; In a 
subsequent encounter this transferred material may become a 
loose wear fragment; / 

Experiments have shown that most of the times the 
softer asperity looses material to the harder one but the 
probability of harder body loosing material is also finite: 
The reason for this is the presence of weak spots within 
the harder asperity as well as its fatigue caused by the 
interactions. In general, therefore, the fracture occurs 
more often on the sicfe of the softer asperity as against 
that of the harder asperity; 
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Brockley and Fleming [13] conducted the experiments 
on copper metallic junctions cut from one piece. These 
single piece junctions were made to shear by applying a load 
from a motor. They studied the failure of these junctions 
in detail. The failure of these junctions is shovm in 
Fig. 1,4. The figure illustrates that a junction formed by 
two identical materials has equal probability of failure on 
either side of the interface. This figure also indicates 
that compression, rotation, tension and shearing at the ends 
of the junctions are apparently involved in the process. 

Thus the mechanism of wear fragment formation is quite 
complex and does not render itsej-f to convenient analytical 
solution. However rough estimates of the forces exerted 
through such junctions have been attempted by Green [14], 

Edward and Hailing [22], and Gupta and Cook [23]; 

New let the Fig;1;6(a), where the asperities have 
just met, be considered. A mean normal stress (compressive) 
and a tangential (shear) stress act through the junction; 
Due to these stresses the junction is defoimed so that the 
two surfaces move parallel to one another, These stresses 
calculated by Green [14] are shown in Fig;1;7. For normal 
values of surface roughness, angles 6 ^ and 0^ shown in 
Fig;1.3 are less than 10°, when a junction is formed. For 
such values of 6^ and it can be seen from Fig. .:1 ;7 
that 1.25 and ^ 2 ICj, v^ere represents the 

yield stress in shear of the softer body. Further deformation 
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as shown in Fig;1:6(c) and Fig;1 ;6( d) , causes changes of 
shape such that increases 6 ^ decreases. It is also 
seen from Fig.l ;? that as decreases and increases 
tends to and decreases continuously. Also as 
the junction tends to become roughly symmetrical continues 
to decrease. Thus, when the junction has been deformed to 
the symmetrical state the norrml stress on it is zero and 
the tangential stress is nearly equal to its initial value. 
After the symmetrical shape has been achieved the cycle gets 
reversed, i.e, increases slightly while P^ becomes tensile 
and remains so until the junction eventually breaks under 
combined shearing and tensile forces. Thus after the 
asperities have matched up and firmly welded the cross section 
’a’ remains approximately constant until necking and fracture 
commences ; 

The formation of a loose wear particle can be under- 
stood to occur as a result of the following stages; 

(a) Matching up of two asperities 

(b) Adhesion of these asperities across the interface 

(c) Fracture initiation at a critical point 

(d) Transfer of material from one asperity to 
another, and 

(e) Subsequent disintegration and removal of ifiis 
transferred material; 
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1 *4 Objective and Scope of Present Work ^ 

In this present work an attempt has been made to 
find the path of fracture at welded asperity junction 
formed during adhesion of two bodies in sliding contact. 

The path has been obtained using finite element analysis 
of the junction; Since FEM is a powerful technique, it 
was considered that it could be applied to fracture of 
asperity junction. The possible advantage could be that 
FEM can easily account for material nonlinearity, different 
geometries of the asperity junction and complex boundary 
conditions ; 

The present work is basically divided in two stages. 
The path of fracture for two different kinds of asperity 
junctions has been studied, (l) For a brittle material 
junction wherein the fracture takes place instantly without 
the junction going into the plastic deformation stage; This 
analysis was done using elastic analysis of finite element 
technique; (2? For a ductile material junction wherein large 
plastic deformations occur before fracture; Here the Elasto- 
Plastic Analysis of Finite Element Technique was used. This 
non-linear analysis was done using the initial stress method; 

Though the present work was basically done to trace 
the fracture path in an welded asperity junction, it can be 
very easily used in any problem where the path of fracture 
is of interest. As an example a notched plate subjected to 
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uniform tensile loading on one edge and fixed along the 
other edge has been analysed to show the path of failure, 
similarly any type of structural problem can be solved using 
the computer code developed* Due to the versatality of 
the finite element technique a welded asperity junction of 
two different materials can also be solved with a few minor 
changes in the codei 

The. main scope of this work would be the application 
to estimate crater wear in cutting tools during machining, 
where the crater wear occurs due to basic adhesive wear; Wear 
studies on bearings can also fall in this category. 

The original idea was to do extensive analysis of 
asperity failure of bimetallic junctions. Due to several 
limitations the work was confined to the following: 

(a) Testing of FEM Model for a brittle material 

« 

junction using the elastic analysis of finite' 
elements. 

( b) Analysis of fracture at the interface of a 
ductile junction using the elasto-plastic analysis: 

(i) when fracture occurs at the interface in 
only one asperity 

( ii) when fracture occurs at the interface in 
both the asperities forming the junction. 
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Qiapter 2 contains a brief description of the 
finite element method. A non-linear FE technique using 
initial stress method is also described; 

Chapter 3 describes the problem formulation and how 

exactly the cieck path is traced making using of the finite 

element method. 

# 

In Chapter 4, the paths obtained using the FE analysis 
are given along with a brief note justifying the acceptability 
of the model proposed; 



CHAPTER 2 


FINITE ELEMENT APPROACH 


2.1 Introduction ; 

The FEM' usefulness and its versatality lies 
in its ability to easily accommodate material nonlinearities, 
geometric non linear! ties and to account for complex boundary 
conditions; This chapter briefly outlines the theory of 
finite element method and its formulation; Displacement 
formulation of the element suitable for the analysis of crack 
propagation in a welded asperity junction is described using 
a 6-node isoparametric element. The initial stress method to 
solve the elasto-plastic problem is presented together with 
a brief discussion of the yield criteria. At the end, a 
brief note on the ccxivergence of the Finite Element Method 
is given, 

2.2 The Finite Element Method ; 

Since FEM is a powerful technique it was considered 
that it could be applied to fracture of asperity junction; 

The possible advantage could be that FEM can easily account 
for material nonlinearity, different geometries of the 
asperity junction and complex boundary conditions; In recent 
years, it has been applied extensively in many areas such as 
structural mechanics* soil mechanics, heat conduction etc. In 
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many situations an adequate engineering model of the system 
is obtained by replacing the system by a finite number of 
well defined components called as elements. Such problems 
are called discrete problems. Even if the number of elements 
is very large, the discrete problem can now be solved readily 
with the advent of digital computers . The various discreti- 
zation methods suggested from time to time to solve the 
realistic Continuum problem, both by Engineers and Mathematicians, 
all involve approximations which approach the true continuum 
solution as the number of discrete variables increase; 

For problems involving complex material properties 
and boundary conditions the engineer resorts to numerical 
methods that provide an approximate but acceptable solution; 

In most numerical methods the solutions yields approximate 
values of the unknown quantities only at a discrete number of 
points in the body. This process of selecting only a certain 
number of discrete points in the body is termed as discreti- 
zation. 

Finite Element Method is now very well described in 
a number of books [24,26] . FEM is a numerical discreti- 
zation procedure by the use. of which a wide range of complex 
boundary value problems can be analysed. Elastic, non-linear 
elastic and elasto-plastic constitutive relations can be 
implemented within the finite element framework in a straight- 
foiward manner; 
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In FEM, the region of interest is divided into a 
finite number of simply connected sub-domains or elements. 

An approximate functional value of the solution is assumed 
over the sub— domain so that the parameters, say Uj^ ( i = 1,i**n) 
of the function becomes the unknovm of the problem. The 
elements are assumed to be interconnected at discrete number 
of points situated on the boundaries and these points are 
known as nodes. The nodal displacements will be the basic 
unknowns, A set of functions is chosen to define uniquely 
the state of displacement within each element in terms of 
nodal displacements; The displacement functions now define 
uniquely the state of displacement within each finite element 
in terms of nodal displacements. Strains along with the 
constitutive properties of the material will define the state 
of stress throughout the element and hence along the boundaries. 
This is discussed in more detail in Section 2;5; 

FEM is endowed with two basic features which account 
for its superiority over other numerical methods. First, a 
geometrically complex domain of the problem is represented 
as a collection of geometrically simple subdomains called 
finite elements. Second, over each finite element the approxi- 
mation functions are derived using the basic idea that any 
continuous function can be represented by a linear combination 
of algebraic polynomials. The approximation functions are 
derived using concepts from interpolation theory and are 
therefore called interpolation functions. Thus the FEM can be 
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interpreted as a piecewise .application of the variational 
methods, in which the approximation functions are algebraic 
polynomials and the undetermined parameters represent the 
values of the solution at a finite number of preselected points 
called nodes; on the boundary and in the interior of the 
e lament ; 


2»3 Basic Approach of FEM ; 

The basic concept is derived fran stoructural analysis; 
Here every structure is approximated as a physical assemblage 
of individual stiectural components or finite elements. The 
elements are interconnected at a finite number of nodes and 
sometimes along the boundaries of the element; Assuming the 
approximate behaviour of individual elanents, the behaviour 
of the entire system can be analysed; After assembling the 
individual elements* the necessary boundary conditions are 
imposed, with primary nodal values as unknowns; Solution of 
the resulting set of equations yield the system response; 


Considering the nature of each individual component in 
an element, the relationship between the primary variable u® 
and the forcing function Q®, for a linear system can be 
written as [24] 




e 


repetitive indices indicating summation; 


( 2 . 1 ) 



27 


The connection among elements can be established 

us ing, 

(a) One set of variables, Uj, for the assembled 
system, i.e, the condition of continuity 




0 

(b) The equilibrium of the second set of variable Q^, 
at each node and equating it to zero. 


NEL 

S 

e=1 



0 


( 2 ; 2 ) 


in which, NEL, is the number of elements considered using 
the above three equations the resulting equation can be 
written as 



[K] 

t u} = 

{ Q } 


(2;3) 

in which 



NEL 




[K > 


2 

e=1 


(2.4) 

and 

1 Q } = 

°oi = 

NEL 

2 

e=1 

CD O 

a 

(215) 


2.4 Finite Element Displacement Formulation of an Elastic. 
Continuum : 

The special features are as given below [ 25] ; 

The continuum is divided by imaginary lines into a 
finite number of elements} these elotnents are assumed to be 
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interconnected at discrete number of nodal points along their 
boundaries. The nodal displacements form the basic unknowns; 

A set of functions are chosen to define uniquely the state of 
displacement within each element. The displacements now define 
uniquely the state of strain within each element, strains 
together with the material properties define the state of stress 
within the element and along the boundary of the element; 

Once this stage is reached, the solution procedure 
follows the standard discrete system pattern; The above 
procedure involves a series of approximations. Ensuring 
displacement compatibility between adjacent elements may not 
be possible always, though within each element it is obviously 
satisfied. Concentrating the equivalent forces at nodes, 
equilibrium conditions are satisfied in the overall sense 
only; 


The displacements can be expressed as 


a. ® f 
1 


u =2N. = [N^, N 


= Na 


( 2 ; 6 ) 


where U = displacement at any point within the element 
N = prescribed function of position called shape 
function 

a® = nodal displacement vector for a particular 
element; 
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The strain vector, e , can be obtained with displace- 
ments known at all points within the element, 

{e} = [ B ] U } (2.7) 


For a linearly elastic constitutive law, 

{a} = [ D]Ue - e^)}+ {a^3 (2.8) 


where, 

[D] 


=s initial strain matrix 

= stress strain matrix 
= initial residual stresses 

nodal force vector which are statically 
equivalent to the boundary stresses and 
.,/ distributed loads on the elements; 



2;5 Displacement Formulation of the Element Used : 

The unique description of the displacement within 
each element in terms of nodal values at boundary points or 
internal points of the elemait is the basic step in any 
displacement finite element formulation and Can be expressed 

clS 

c u } = [N ]Ca®} (2;9) 

where [Nj is the matrix of shape functions; 
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If the shape function chosen to describe the element 
geometry are identical to those used to prescribe function 
(displacement) variation, then the element is teamed isopara- 
metric. The basic procedure in the isoparametric foomiulation 
is to express the element coordinates and element displace- 
ments in the foimi of interpolation functions using natural 
coordinate system of the element. The functional property 
of the interpolation function is that its value in the 
natural coordinate system is unity at node i and zero at 
all the other nodes. The formulation given below follows 
the general pattern of derivation suggested by Bathe and 
Wilson \26l and Zienkiewicz [24] ; For a 2-dimensional 
toriangular element , interpolation functions for a 6— noded case 
is given below, can be derived from the general quadratic 
relationship 


Ni 


h + 4^’ Lg + L 3 + 


+ a^^a^L. 


5 ^ 2*^3 


qI i) L L 
^ H ^3 

( REFER APPENDIX aVi ) 


( 2 : 10 ) 


A quadratic variation of displacement and geometry 
are assumed for this elemaat and it has 12-displacement 
degrees of fooeedom (Fig; 2:1)* 

Similarly, the displacements expressed in teonns of 


nodal values are: 
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u = N,u, + + N 3 U 3 + (2;11) 

V = N,v^ + NjV^ + N 3 V 3 + + N 5 V 5 + (2'i12) 

along the edge of theelement the variation of displacement 
is quadratic* 

Writing equations 2,11 and 2,12 in matrix notation 
tu1 = [N] Ca } (2:13) 

where 

(aJT =[u^ V, U3V3 U4V4 U3V3 u^Vg] . 

The strains at any point within the element is obtained 
knowing the nodal displacement, by, 

ie} = [B]{a} (2:14) 

where [B] is the strain displacen^nt matrix (REFER APPENDIX A:2) ; 

The stress at any point within the element is now 
obtained knowing the strain vector and the stress strain 
matrix [D ] by 

ia} = [ D] {e} (2.15) 

where [D] is the stress-strain matrix for plane stress 
(REFER APPENDIX A:3): 

The matrix to be evaluated are the stiffness matrix 
[KJ;- [IC]is evaluated by numerical integration of 
t / B^D B dA (2:16) 
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/ 



n 

I J I t 2 

i=1 


w. 

X 


f (L’ 


CD 




(2;17) 


n stands for the number of Gauss point. 


matrix. 



9 x/dL^ 9 y/5L^ 

3 x/^L 2 9 y/ai^ 


J is the Jacobian 

(2;18) 


weighting functions of gauss point where n for 
cfuadratic element is 3. 2ULenkiewicz [24] has given the 
following weighting functions for quadratic triangular 
element using Gaussian quadrature formula with 3 points; 


2,6 E las to- pi as tic Analysis ; 

The difference between elastic and plastic behaviour 
under uniaxial sttess is brought out in Fig;2.2. In non- 
linear elastic behaviour the stress can be expressed as a 
function of strain as a =o’ (e). The main difference of 
plasticity formulation from nonlinear elastic formulation is 
such an explicit relationship is not available; Althou^ 
the stresses at any level of strain have to lie on or within 
the current yield surface, the exact value of each component 
cannot be determined; 

Many numerical solutions of elasto-pla stic problems 
by FEM have been developed, the most common are as listed 
belcw; 

(1) Incremental tangential stiffness approach} 

(2) Initial strain method} and 

(3) Initial stress method [27 ] ; 
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a 



(j 



c. Strain hardening 
plasticity 


a 



a 



d. Ideal plasticity 


Fig. 2.2 Elastic and plastic behaviour under 
uniaxial stress 
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Initial strain method fails to converge for perfect 
plasticity and a small degree of work hardening because of 
the large plastic strains that occur. Initial stress method 
overcomes this difficulty because it relies on the fact that 
a unique stress exists for a increment of strain; Hence 
the last one. is more popular and is discussed in detail in 
this chapter; 

2.7 Yield Criteria ; 

A yield criteria is a hypothesis concerning the 
limit of elasticity under any possible combination of stresses. 
The suitability of any proposed yield criteria must be checked 
by conducting experiments. It has been shown that the yield 
'criteria is a function of - ^ 2 )* ^^2 ”^ 3 ^ 
and is independent of the hydrostatic stress ccxnponent 

’ ( +^2 

Thus yielding occurs when some scalar function of 
the principal stress differences reaches a critical magnitude, 
mathematically: 

f(ai - cr^ , 0 -^ = constant (2;19) 

f 

The two of the most commonly used yield criterias are: 

( 1 ) Tresca Yield Criteria 

(2) Von-r-Mises Yield Criteria; 

Von-Mises yield criteria used in this work is 


discussed below: 
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Von-Mises Yield Criteria : 

Mathematically it can be stated as, 

^ ^ 3 ” ~ constant 

( 2 ^: 20 ) 

In this function, each of the principal stresses 
contribute to yielding; The criteria was interpreted by 
Hencky as "Yielding begins when the shear strain energy 
reached a critical value" ; 

The value of the constant in the above criteria can 
be determined from the simple tension yielding^ i;e. = ^yp^ 
^2 ~ ^3 ~ ^ pure shear yielding, i;e. “■'®3 “ 

cr^ *= 0. Hence the Yield Criteria becomes. 






2 2 
2CT = 6K^ 

yp 


( 2 ; 21 ) 


where cr ^ “ tensile yield stress and K is yield shear stress, 

y H 

For a general state of stress i ^ -^t o o . cr , 


^ ) , The Von-Mises yield criteria can be expressed as 


<N 

+ 


CM 

+ c®. 

^ + T 

2 

+ ’ 

ni/2 

= cr 

xy 

yz 

ZX' 

'j 

i 


( 2 : 22 ) 


Figure 2,3 shows both the yield criteria represented 
in principal stress space and TT plane, respectively; 




P’s 



(b) 

Fig. 2.3 Stress space representation of yield criteria 
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ElQS'bo— pXb S' ti.c S'tresg—s'trctin Matrix ? 

Hookes Law for the isotropic elastic material can 
be written in the matrix form as ia}- [DjCek 

Yaraada and Yoshimura [28^ have developed a plastic 
stress-strain matrix [Dp^J (SEE APPENDIX A:4) by inverting 
the Prandtl-Reuss equations obeying the Von-Mises yield 
criterion; 


For the plane stress case the [D ] matrix will 
reduce to the following as given by Yamada and Yoshmura [28] 

‘2 

V ^ 


D = 
o 


where j 


E 

Q 


(- \ + 2yp) ( + 2p) 


(- 


cr' + vcr' r ) 
X y xv^ 


1 +V 


(-a +2^0* T ) 

^ y X xy^ 


Symmetric 


9E 


2H^ ( 1-v) 


o - ^52 

P - +1-+» 




(2;23) 


Q = R + 2 ( 1 -V ) p 


10 • < 10 

R = cT + 2^cr a +0 
X X y y 


( 2:24) 
(2;25) 


cr„ and stand for the deviatoric stresses; 
The elasto-plastic matrix is then given by 


[Dep ] = [D] - [ Dp] : 


( 2 I 25 ) 
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2S The Initial Stress Computational Method : 

In the initial stress method [ 27] the solution of 
the non-linear problem is approached in a series of approxi- 
mations. During a load/ displacement Increment, a purely 
elastic prdclem is solved, determining an increment of 
strain {Ae‘} and stress {Act'} at every point in the 
continuum; For the increment of strain found from the analysis, 
the stress increment in general will not be correct. If ^ 
is the true increment of stress possible for the given strain 
then the situation can only be maintained by a set of body 
forces, equilibrating the initial stress system given by 
{A o 1 - {Act} . At the second stage of the computation, this 
body force system can be removed by allowing the structure 
(with unchanged elastic properties), to deform further. An 
additional set of strain and corresponding stress increments 
are caused. Once again these are likely to exceed those 
permissible by the non-linear relationship and the redistri- 
bution of equilibrating body forces has to be repeated. If 
the process converges then within an increment the full 
non-linear, compatibility and equilibrium, conditions will be 
satisfied. In each cycle, since the same elastic problem 
is solved only partial inversion of the elastic equations is 
needed; 

For the elastic-plastic situation, the steps during 
a typical load/ displacement increment can be summarised 
as follows [27] t 
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1, Increment of load/displacement is applied. Elastic 

increments of strains fAe*} ^ and corresponding elastic 
increments of stress are determined^ 

{a cr* is added to stresses existing at sta-rt of 
increment, £ }, to obtain {o’}, and thus *ar . Check 

whether a < 'the condition is satisfied 

only elastic strain changes occur and the process is 
stopped. If the condition is not satisfied proceed 
to the next step; 

— i. * 

3. If o’ and if (the element has 

yiel^d at the start of iteration), £ Acr}^ is found 
by £Ao}^ =[^Epl (2;27) 

where Dgp is computed using {cr*}; Stresses which 
are supported by body forces are evaluated 

-r 

{Aoll “ £Ar*} -j ''£Acr}| (2;28) 

Current stress is 

to} = { cr •} - CAer"} ^ (2;29) 

Current strain is 

{el = £e '} + £ae5^ (2.30) 

W 

4, If o ’ > Oyp, but i;e;, elastic-plastic 

transition has occurred in that iteration, a factor 
H from equation (2;31) is found and then the 
intermediate stress value is found out at which 
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yielding begins. The increment in stress is computed 
by the relation, { =: [^EP^ from that 

point. Then proceed as in step 3 with computed 

from stresses 

_ 1 

Factor = — I — (2;31) 

a “ cTq 

5. Nodal forces are computed corresponding to the 

equilibrating body forces. These are given by 

(RlJ® =^^[6^] {Aa^ d (vol) (2;32) 

elemental nodal force are added to get the complete 
nodal force vector tR]|} : 

6; The problem is resolved using the original elastic 

properties and the load system to find fAo’’}^ 

and { A e '}2 • 

7; Current values are found., 

8. Steps 2 to 6 are repeated; 

The cycling can be stopped when the appropriate 
convergence criterion specified is satisfied.. 

2.10 Convergence of the Iterative Procedure ; 

The adoption of an incremental plasticity theory 
and iterative solution technique, such as used in "Initial 
Stress Method" leads to errors due to size of increment and 
finite values of criteria specified in a solution, indicate 
whether the non-linear equations describing the behaviour of 
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the material when subjected to increments of load/ displace- 
ment have been solved with sufficient accuracy. This is 
achieved by comparing the change in some aspect of the 
solution between iterations and terminating the procedure 
once this change has become sufficiently smalli 


Nayak and Zienkiewicz [29] have described some of 
the convergence criteria for load based analysis and displace- 
ment based analysis. The convergence criteria imposed in 
this work is on the change in the absolute magnitude of 
largest term in load correction vector between the iterations, 
compared as a percentage of the absolute magnitude of largest 
term in load correction vector at that instant calculated 
by equation (2.33), 



(2.33) 


where, 

C =: convergence criterion 

= magnitude of largest terra in residual load 

nmax i 

vector at ith iteration 

Rj^ = magnitude of largest term in residual load 

max i +1 

vector at (i+1) iteration; 

The following chapter briefly outlines the problem 
formulation and how the path of crack is determined using 
FEM. A note gri formation of double nodes as crack propagates 
is also given; 



CHAPTER 3 


PROBLEM FQRMLJLATIQN AMD METHODOLOGY 


3,1 Introduction : 

This chapter describes the details of the actual 
problem of welded asperity junction and the finite element 
methodology of tracing the path of crack in the asperity 
junction. The asperity junction fracturing is considered 
as occurring under plane stress conditions: Six-node 
isoparametric triangular elements are used for the analysis 
and numerical integration is employed for calculating stresses 
and strains; A list of the assumptions and the description 
of the data vtiich is necessary for the execution of the 
programme is also mentioned, together with the solution 
procedure. The criteria used to determine crack initiation 
and propagation is also stated. 


3,2 Problem Description : 

The asperity junction is considered at a stage of 
its life cycle wherein it is symmetrical as shown in Fig; 3.1 
For the present analysis, the function as it fractures is 
simplified as occurring under plane stress conditions; In 
this work an attempt has been made to obtain the distribution 
of stress and strain at the nodes which ia then used to trace 
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FIG.3!|I)FINITE element mesh used 


^ 4-6 
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the path of crack; 

The process of junction fracture is analysed using 
the finite element method with the following assumptions: 

( i) The junction geometry is assumed to be 
perfectly symmetrical about the interface. 

( ii) The junction failure is assumed to occur under 
conditions of plane stress. 

( iii) Small deformation theory is applicable for 

the purpose of calculating stresses and strains. 

( iv) For the ductile junction considered the material 
considered is isotropic, linearly elastic, 
linear strain hardening type obeying the Von- 
Mises yield criterion. 

The following geometric boundary conditions and force 
conditions are specified for the present work as shown in 
Fig; 3.1 ; 

(i) The boundary AB is fixed both in the X and Y 
directions. 

(ii) The boundary EF is subjected to force boundary 
conditions in the X direction; 

3,3 FE Formulation and Crack Propagation : 

The problem described above is solved using FEM 
explained in Chapter 2, The asperity junction is discretized 
as shown in Fig;3;l . The study of wear, involved the modelling 
of crack propagation usually emanating from the asperity junction. 
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The propagation is simulated by using the concept of double 
nodes. 

The crack propagation in finite element method is 
simulated by making each node along the crack path into two 
nodes (double nodes) once the crack passes through that node; 
Consider a path along which crack is going to propagate as 
shown in Fig;3;2. As the crack passes through the node i, 
it has to be replaced by two nodes and as the crack passes 
through the node j another two nodes have to replace the 
node j, because the elements I and II are now no longer 
connected at i S<j as was the case before the crack propagated 
along that edge. Due to this formation of double nodes 
the node numbering has to be redone and a brief description 
on how it is done is given below i 

The elements connected to node i (before crack 
propagates) are all identified and the slopes made by node i 
joining their (elements) midnodes on edge opposite to i 
are found out and stored. Now the angle made by the crack 
along the edge of element to reach the node is found out, say , 

this angle is with reference to node i; The path along which 
crack now is going to pass is identified and the angle made 
by that path with respect to node i is again identified, say 
At this stage the node number of i and j are compared as 
to find out which one is greater. With this information 
the node renumbering can be done as follows! 
if i > j, then 
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( i) for all elements connected to node i and 
slope made by mid point with i; between 8^ 
and ©2 increase i to i+1 • 

( ii) for all elements connected to node i and slope 
made by mid point with i; not between 0^ and 
©2 increase i to i+2; 

( iii) for all elements having node numbers greater 
than i increase node number by 2: 

( iv) for all nodes having node number between i and j 
increase node number by 1 ; 

(v) for element containing node j if slope made by 
mid point with i between and © 2 f no change 
in nodfe:n umber ji 

(vi) for element containing node j and if slope not 
betottom©^ and ©2 increase node number to j+1 ; 

(vii) for nodes having number below j no changes; 

Similarly a set of rules for renumbering, if j > i ; 

Because of the crack the nodal numbers change and 
due to this a new connectivity of elements has to be 
generated to account for the crack using the node numbers 
as given by the above set of rules. However it is to be 
noted here that the coordinates of the double nodes are 
assumed to have the same coordinates as that of the parent 
node. Thus it is clear that though the node nunfcer and 
connectivities have changed the elemental stiffness matrices 
stored on tape are unchanged. Using these elemental matrices 
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and new connectivities generated the global stiffness matrix 
is reassembled* 

Thus it is seen that crack propagation requires 
( i) double node formation; ( ii) node renumbering and 
(iii) corresponding change in connectivity and ( iv) reassembly 
of global matrix; 

All this has been incorporated in the computer code 
developed: 

3*4 Direction of Crack Propagation ; 

From the finite element programme the forces, displace 
ments, stresses and strains within the junction are found out. 

From the experimental studies it was observed that the crack 
initiates at the interface, using this, (i.e, knowing the 
particular node at which the crack initiates) all the paths 
(i.e, sides of elements) along which the crack can possibly 
propagate are investigated to see if the criteria for crack 
to occur has been satisfied; 

The procedure for checking the criterion is briefly 
outlined below. Consider a small section (Fig*3, 3) wherein 
the crack has propogated till node ”B*', Now the possible 
paths for the crack to propogate are BC, BD, BE, BF, BG, BH, 
and the possible elonents are I, II, III, IV, V, VI, VII, 

First of all the elements along wbich. the crack can possibly 
occur are identified in the continuum, along witii the identification 
of the elements. The possible paths and the coordinates of the 
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corresponding nodes involved (i,,ei end nodes of path and 
mid node along path) are found out. The position of the 
elements around the node "B" are identified by finding 
out the slope made by line joining "B" to raid node of the 
side opposite B (eg. line B1); Also, the slopes made by 
the possible paths are also found out and s.tored., , Now, if 
the crack was to propagate along BE then the shear stress 
along the edge BE of Element III would have to exceed the 
maximum shear stress of the material; To 

find this tangential stress, firstly the force contribution 
at node B from all elements lying between are summed up, 
the component of the force along the line BE of this force is 
found out and this is then divided by the length of the 
side BE to give the tangential stress along ”BE", If the 
crack criterion is satisfied then the crack is assumed to 
have proceeded along BE. Now, the node number along the 
path of crack are noted and using these the double nodes 
formed due to the crack are correctly numbered and all 
previous node numbers above the cracked node number also 
undergo changes. Because of this change in node numbers 
the global stiffness matrix at each stage will be different. 
Hence, using the stored elemental matrix we can again assemble 
the global stiffness matrix in accordance to the new connectivity 
generated within the programme; From this step the process 
is again repeated as above till the other edge of the junction 
is reached; 
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3 -5a Methodology for Brittle Materials: 

Using the FEM the stresses at all points in the 
continuum are got and knowing the node from which the crack 
will emanate all possible paths from that node are identified 
and checked to see if crack can propagate along it as mentioned 
in Section 3.4. Node renumbering etc. as mentioned in Section 3,3 
due to double node formation because of crack is done. Again 
the whole process is repeated till the other edge of the 
junction is reached; 

3;5b Methodology for Ductile Junction ; 

The procedure for the analysis of crack path is 
more or less similar except for the fact that here nonlinear 
FEM analysis using the "Initial stress" method is used. 

Here, the Icfad is applied in steps, and at each step the 
possible crack paths are analyzed and checked if crack criteria 
is satisfied. If the crack has propagated the procedure as 
mentioned fbr brittle junctions, i;e. from changing nodal 
numbers to again checking for crack propagation is repeated. 

In case the crack does not propagate along any path then the 
next load increment is applied and so on. 

In this incremental cum iterative technique of elastic- 
plastic finite element analysis the size of the increment 
and the number of iterations to be performed within an 
increment will affect the results obtained. For this work, the 
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values of load specified ij\ the first load step is such that 
yielding just begins. It has been observed that by using , 
the initial stress method for the nonlinear material problems 
even for large values of load/ displacement increments, the 
results obtained are quite satisfactory. Because of the 
large computer storage required and large execution time of 
programme the number of load steps has been kept very small 
in the present work, 

3 .6 Input Data Structure : 

The following data was supplied as input for the 
purpose of analysis: 

(1) Total number of nodes and elements, weighting 
function for numerical integration, maximum number of 
iterations for each load increment, bandwidth of matrix, etc« 

(2) Material properties like values of elastic 
modulus, Poissons ratio, yield stress, ultimate shear stress, 
strain hardening coefficient; 

(3) Nodal data this includes the node number and 
the x,y coordinates of the node for all the nodes. 

(4) The element data in which the element number is 
given and the nodal connectivities of the element, for all 
the elements. 

(5) The boundary node number on which the boundary 
condition is known and the type of boundary condition on that node; 

A brief description of the results obtained and the 
correlation with some known results are briefly described in 



CHAPTER 4 


RESULTS AND DISCUSSICa^IS 


4,1 introduction ; 

This chapter presents an analysis of the results 
obtained for the fracture path of the welded asperity junctions; 
The finite element analysis has been done using three type of 
meshes. The coarse mesh had 213 nodes with 88 elements and 
the fine mesh had 253 nodes with 108 elements. The CPU time 
required for one run of the elasto-plastic analysis on an 
average has been around 12-15 minutes on the DEC 1090 computer 
system; 

The elasto-plastic computer code developed has been 
compared with a cantilever beam problem [27] • Load versus 
displacement diagram has been plotted (Fig;4;1) and the 
results compared fairly well with that of [27] ; Zienkiewicz [ 2 ^ 
has used 150 triangular elements with a total of 99 nodes, 
whereas the present analysis used to check the computer code 
had 24 six-node isoparametric elements with a total of 63 
nodes; 

Figure 3;1 shows the finite element discretization of 
the asperity junction together with the geometric bouncfery 
conditions applied. The meshes have been made manually and 
finer meshes were not tried out, because that would have 
involved a higher CPU time; 
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4.2 Brittle Junction Failnrfi ; 

In the case of asperity junctions of the two brittle 
materials as mentioned earlier, the elastic analysis of the 
continuum was done. The calculations showed that fracture 
started at edge A (Refer Fig; 4;2) on the interface and 
propagated along the path A-1 , 1-2, 2-3, 3-4: . At 

each stage crack has been assumed to propagate which the shear 
stress along the edge of an elonent exceeds the ultimate shear 
strength of the material; 

Thus the ccmputer results obtained showed that the 
fracture path, which initiates at the junction edge A, has taken 
a straight path along the interface. This result is in agreement 
with the intuitive feeling with which the brittle junction 
case was used to test the FEM model; 

The intuitive feeling was also supported by the obser- 
vation that in case of brittle materials like cast iron rubbing 
against each other there is very less wear compared to the 
case of steel and other ductile materials which have good 
adhesive properties. If the fracture of welded junction in 
brittle materials would be occurring in a similar manner to 
that of ductile material the (adhesion) wear coefficient would 
not have been drastically lower for similar conditions. But 
several published results [15,30] have shown that wear of brittle 
materials is indeed much less; 

The adhesive wear is shown to be governed by the equation 

P 

u 


V 


2 


( 4 . 1 ) 



mm 
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where V 
L 
Z 
P 
H 


Volume of wear 
Distance of sliding 
Wear coefficient 
Contact load 

Hardness of the material ( softer) ; 


when the two bodies are under contact the load is supported 
by some asperities which ultimately form junctions; The 
number of junctions formed can be estimated by an approximate 
expression [15 j 


P = HN a (4;2) 

where a = Area of contact at each asperity junction 
supporting the load 

N = Total number of asperity junctions supporting 
the load; 

Table 4;1 shows values of Z for four different materials 


TABLE 4.1 [30] : 


Material 

— 

Z X 10“^ 

Remarks 

I . - - 

Hardness 

Zinc/ zinc 

160 

High adhesion 

Low 

Cu/Cu 

32 ’ 

Medium adhesion 

Medium 

MS/MS 

.01 

Lew adhesion 

Hi^ 

we/wc 

.001 

Very low 
adhesion 

Very hi^ 
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Referring to Fig;4,3 it can be observed how the 
wear rate increases as the number of asperities (N) 
increases [31] . It can be also seen that as the hardness 
increases the number of asperities decrease and the wear 
rate also decreases; Now, since in engineering materials 
hardness and brittleness have a direct correlation, it is 
obvious that adhesion wear of hard and brittle materials 
is mainly by abrasion and not by adhesion. The adhesion 
characteristics of hard and brittle materials indeed 
decreases; Thus, it may not be wrong to assume that fracture 
in a brittle juncticn may be straight one as A-1 7- B ; 
This qualitative comparison then lends support to the 
inference that fracture path A-1, 1-2, 2-3, 3-4, 4-5, 5-6, 
6-7, 7-8, 0Bi refer Fig, 4.4 is correct and hence was 
considered acceptable for further application; 

4,3 Ductile Junction Failure; 

Trial computations were done to arrive at an approxi- 
mate ultimate load for the asperity junction (Poq)* 

For a ductile junction the elasto-plastic analysis 
using initial stress method was carried out and it was 
found that with fine mesh the fracture initiated at a load 
of 0;6 At 0;6 the crack propagated from A to 1 

as shown in Fig, 4,4; After making the necessary changes 
due to double node formation it was found that the crack 
did not propagate and hence another load increment of 0.1 P 
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was applied and the check for crack propagation was 

applied. It was found after successive load increments 

that at 0;8 P_„ crack propagated along one element 1-2 
ap 

and on checking for further propagation it was seen that 

2-3 also cracked. Finally at load the crack propa- 

ap 

gated from 3 to B along the paths shown; This type of 
crack propagation is expected in the case of bimetallic 
junction; When sliding occurs at a bimetallic junction 
the failure cannot occur in both the bodies simultaneously; 

It will occur in the weaker of the two bodies; This is 
a physical observation true for all rubbing bodies* 

It was observed that after point 4 (Fig; 4;4) 
the crack path showed a tendency to fall to the interface. 
This could be because of the symmetry of the junction being 
lost due to the opening of the crack, A plot of load 
applied versus crack length is given in Fig. 4,5; 

Crack propagation in the coarse mesh was also 
analysed and the results are shown in Fig; 4,6; When a 
junction is formed between two asperities of same materials, 
the fracture can occur with equal ease in both bodies. 

Thus the initiation of the fracture must occur simultane- 
ously as already discussed in Section 1.3 , This indeed 

has been verified by Brockley [13] through his experimental 
(Figure 1,4). Thus having checked the finite element model 
for a bimetallic junction and a brittle material-j unction 
the model was applied to analyse a junction of same material; 




FIG 4 . 6 . DUCTILE JUNCTIONT FAILURE COARSE MESH 
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Therefore, fracture path of asperity junction with crack 
starting from both sides was carried out by incorporating 
a few changes in the computer code. Here the fracture 
started in both sides at 0;6 Again, after first 

element edge had cracked the crack did not propagate and 
only at 0.8 P crack started propagating from both 

Cl 

sides; But it was observed that after crack had propa- 
gated along two more element edges it was found that the 
rate of convergence became very slow and large number of 
iterations would have been required before the criterion 
for convergence was satisfied. Hence the analysis was 
stopped at this stagei 

The initiation of cracks from both sides agreed well 
with the experimental work of Brockley and Fleming [13^ 
where their Figure 1 .4 shews that crack follows a path 
very similar to that as obtained using this analysis 
(Fig. 4.7). This, thus supported the acceptability of 
the model suggested in this work; 

Tables 4.2a and 4.2b lend support to the obtained 
results as it can be clearly seen that for bimetallic 
junction the wear coefficient Z is considerably less than 
that for similar metal junction. Thus meaning lesser 
adhesive wear at junction; 




FI6.4-.7 UNIMETTALIG JUNCTION FAILURE 
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TABLE 4.2a [15] Under Identical Conditions 


» 

I Z for metal on metal 

Conditions I 1 h - j - 4!/,V T 4 T 


Different 


Dry surfaces 

5xt0~^ 

10- 

Boundary 

2x1 0""^ 

10- 

lubrication 



Medium 

lubrication 

10-5 

10 -' 


TABLE 4;2b [15] Under Identical Conditions 



I 


Material I Wear coefficient Z 



MS/MS 7 X 10“^ 

60/40 Brass/MS .6 x 10""^ 

70/30 Brass/MS 1 ’;7 x 10“^ 


An approximation to the volxome of material lost 
can be made very easily knowing the fracture path as 
obtained from this analysis; 
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The broad work is however, far frcan complete and 
the model may need to be examined in some detail as 
outlined in the conclusions chapter;. 



CHAPTER 5 


QONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

Cone lus ions t 

Based on the results obtained for the proposed 
model, the following conclusion are drawn; 

(i) FE method can be applied to predict the 
path of fracture' at welded asperity 
junctions ; 

( ii) The maximim shear stress criteria is 
acceptable for predicting thepath of 
fracture* 

(iii) The present approach explains fairly well 
the qualitative behaviour of the fracture 
of asperity junctions in the case of brittle 
and ductile materials i 

The following are a few of the s ingestions which 
the author feels would be of importance towards the 
acceptability of the proposed model, as a standard one; 

( i) Finer meshes should be tried out to fit 
the exact path of fracture in the asperity 
junction; 

(ii) Other criterias like energy release rate, 
stress intensity factors etc; should be 
investigated to find fractusre path; 
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(iii) The author strongly feels that experimental 
work using simulated junction models of 
both brittle and ductile materials should 
be carried out and some statistical results 
obtained* 

(iv) The application of magnetic field on the 
fracture path should also be investigated 
since the influence of magnetic field can 
wear is fairly well established and the 
basic mechanisms like change in fracture 
path is responsible for change in wear 
characteristics; 

(v) The present analysis allows the crack to 
propagate along the edges of the elements 
only*. More realistic analysis can be 
modelled by modifying the analysis so 
that the crack can propagate across the 
elements also; 
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APPENDIX 


N^ = ( 2 L^ - 1 ) 

( 21 ^ - 1 ) 

N3 = L3 (2L3 - 1) (REFER [24]) 

N4 = 4 . 

% = 4 L 2 L 3 

^6 = 4 HL 3 


Cartesian co-ordinates x,y and natural co-ordinates 
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(REFER [24] ) 
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(REFER [24]) 
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(REFER [28]) 


E = Elastic modulus 


, I - , 2 
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^ = Poisson’s Ratio 

S = 2 a-2(i ^hI) 


0 = 


3 

1 


ZX •' 


z' ■ ' z X 

G = Shear modulus 
H' = Strain hardening parameter! 


